Off-lattice noise reduction and the ultimate scaling of DLA in two dimensions 
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OflF-lattice DLA clusters grown with diflFerent levels of noise reduction are found to be consistent 
with a simple fractal fixed point. Cluster shapes and their ensemble variation exhibit a dominant 
slowest correction to scaling, and this also accounts for the apparent "multiscaling" in the DLA mass 
distribution. We interpret the correction to scaling in terms of renormalized noise. The limiting 
value of this variable is strikingly small and is dominated by fluctuations in cluster shape. Earlier 
claims of anomalous scaling in DLA were misled by the slow approach to this small fixed point value. 

PACS numbers: 64.60.Ak, 61.43.Hv 
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I. INTRODUCTION 

Since its introduction in 1981, the Diffusion-Limited 
Aggregation model of Witten and Sander has been 
a paradigm of self-organised scaling behaviour in irre- 
versible growth. However, even after twenty years, there 
is still controversy about its scaling properties; many 
authors have claimed, for example, that DLA clusters 
do not scale as simple fractals, but instead have various 
anomalous features. In this paper we give data on DLA 
clusters with noise reduction which enables us to refute 
conclusively the basis of these claims of anomalous scal- 
ing. We will show that the apparent anomalies arise from 
a slowly decaying correction to scaling which can be as- 
sociated with the level of intrinsic growth fluctuations, 
as suggested in [Q. The analysis of these corrections to 
scaling gives us considerable insight into the asymptotic 
behaviour of DLA, i.e. the DLA fixed point. 

In (off-lattice) DLA a cluster is rigid and stationary, 
growing from one seed particle by accretion at first con- 
tact of N mobile diffusing hard sphere particles. The 
diffusing particles are sufficiently dilute so that they can 
be taken to arrive one at a time. We consider the distri- 
bution of where growth (by deposition) occurs at a given 
cluster size. The average radius of deposition is defined 
by i?dcp = {^), where r is the the distance of deposition 
from the center of the cluster. There is no controversy 
that i?dcp c>c N^/^ , consistent with a simple fractal of 
dimension D = 1.71 for large clusters in two dimensions. 
However the spread of the deposition radius is thought 
to show anomalies. Plischke and Racz Q introduced the 
penetration depth, ^, the standard deviation of radius of 
deposition of a given cluster, and claimed that it scaled 
differently from i?dcp- More recently, Davidovitch et al. 
Q considered the standard deviation of the cluster av- 
erage radius across the ensemble of clusters, 6Ros, and 
claimed that it was asymptotically negligible compared 
to the mean. Another anomalous feature that has been 
claimed of DLA is multiscaling [^|j6j: the fractal dimen- 
sion of the cluster is said to depend on the distance (rel- 
ative to the cluster radius) from the center. We will ex- 
amine these claims using finite size scaling with the help 



FIG. 1. The noise-reduced DLA algorithm. A particle is al- 
lowed to diffuse freely until it contacts a cluster particle. The 
diffusing particle is moved onto the cluster particle, reducing 
the distance between their centers by a factor A. 



of noise reduction and show that none of them hold. We 
find that DLA is consistent with simple scaling, and the 
apparently anomalous scaling can all be explained by a 
slow correction to scaling. 



II. OFF-LATTICE NOISE REDUCTION 

Noise reduction for the lattice version of DLA has been 
introduced |^J^ with the aim of suppressing the shot 
noise of the individual incoming particles. When growing 
at lower noise levels, the clusters achieve more asymptotic 
behavior at smaller sizes: a prime example of this is that 
the lattice effects show up earlier. These lattice effects 
on noise reduced clusters (or without noise reduction on 
very large clusters) are quite strong, so in order to avoid 
them, any analysis of large scale DLA clusters has to be 
made off-lattice. 

In our version of noise-reduced off-lattice DLA the par- 
ticles diffuse freely until they contact a particle in the 
cluster, just as in the original model. However, on con- 
tact with a particle of the existing cluster, the diffusing 
particle is moved into that particle by a factor of A, at 
which point it is irreversibly stuck (see Fig. |l|). This 
means that shallow bumps are added to the cluster, and 
that we must add 1 /A particles on top of one another to 
protrude the growth by a particle diameter. A cluster 
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grown with this method of noise reduction is shown in 
Fig. Another way to do noise-reduction of this type 
was introduced by Stepanov and Levitov , who gener- 
ahzed the method of iterated conformal maps [[lOp|] to 
add shallow bumps. 



III. FINITE SIZE SCALING 

Growing clusters at a variety of levels of noise re- 
duction gives us a very clear picture of the finite size 
scaling effects in DLA. We grew 1000 DLA clusters to 
1,000,000 particles with noise reduction levels oi A = 
0.3,0.1,0.03,0.01, and 4000 clusters with A = I as well 
as 25 clusters with A = 0.001. At various points in the 
growth, 100,000 probe particles were fired at each cluster 
to measure its properties. In the following measurements 
the center of the cluster was taken naturally as the center 
of mass of these probe particles ( "center of charge" ) . 

Fig. ^ shows a primary test of scaling: how the rela- 
tive penetration depth, the ratio of penetration depth to 
mean radius of deposition S = ^/i?dopj varies with N. 
The different levels of noise reduction are all consistent 
with a universal asymptote, Sqo ~ 0.12, and with f^-'^-^^ 
as the common correction to scaling at large N. Fig. ^ 
also shows data obtained with the HL iterated conformal 
map method 0. We can make a naive geometric argu- 
ment to see how the HL bumps correspond to different 
levels of off-lattice noise reduction. In the HL method the 
bumps are generated by a conformal map parameterized 
by a (see |10 for details), with small a giving shallow 
bumps. Using the scheme of Fig. |l|, and working out 
when the aspect ratios of the two types of bumps match, 
we find: (1/a— 1)^ ~ 2/A— 1. Thus a = 0.5 corresponds 
to ordinary DLA. Other equivalent cases are indicated in 
the legend of Fig. |^, showing that this naive argument 
represents reasonably well the relationship between our 
noise reduction and that of Ref . ||^ . 

The correction to scaling exponents we report in this 
paper are not arbitrary fits, but directly measured as 
follows. If we posit a leading asymptotic form of some 
quantity, Q: 



QiN) = Qooil + CN-'') 



(1) 



then a plot of dQ{N)/d\n{N) vs Q{N) should have an 
intercept on the Q axis of Qoo approached with slope 
—V, both independent of the magnitude of C. Fig. ^ 
shows this analysis applied to S and this is the basis 
for the choice of exponent ly = 0.33 for Fig. ||. This 
provides unbiased evidence that all the different levels of 
noise reduction approach the same asymptotic value Sqo , 
consistent with a common correction to scaling exponent. 
Interestingly the 'fixed point' can be approached from 
either side, corresponding to opposite signs of C. 

A deeper test of the universality of these clusters and 
their scaling comes from the niultipolc moments of the 
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FIG. 2. Behaviour of the relative penetration depth 
H = (,/Rdcp, with varying cluster size at various levels of 
noise reduction. The abscissa is chosen accordine to the cor- 
rection to scaling exponent measured from Fig. H. The rela- 
tive penetration depth clearly converges to a non-zero com- 
mon value. Also shown are curves for clusters grown by the 
Hastings-Levitov (HL) method with expected equivalence in- 
dicated in the legend. The top right panel is a magnification 
of the asymptotic end of the curves. 
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FIG. 3. Rate of change of the relative penetration depth 
d'E/d\n{N) plotted against H. The common dashed asymp- 
tote indicates that H has a dominant correction to scaling of 
the form H = E^{1 + CN'") with Hoo = 0.121 ± 0.003 from 
the intercept of the plots and v = 0.33 ± 0.06 from the slope. 
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FIG. 4. The rate of change of the multipole powers P2 and P5 with respect to ln(A'^) against the multipole powers. Each 
show a correction to scaling exponent of 0.33 (dotted hnes), within statistical error. 



growth probability distribution. The nth multipole mo- 
ment is given by 



M,, = dq{x + iyf 



(2) 



where q is the probability distribution for where growth 
will next occur. (Note that q is equivalent to the charge 
density on the cluster surface when it is considered to 
be a conductor held at a fixed potential). The multi- 
pole moments for positive n fully characterise the cluster 
shape, and can be related invertibly to the Laurent coef- 
ficients of its conformal map from the unit circle In 
practice we measured the Mn by sampling (x + ij/)" with 
non-growing probe particles. 

Fig. ^ shows the correction-to-scaling analysis of the 
corresponding multipole powers 



E?2n 



(3) 



where we have scaled each Af„ by the appropriate power 
of the effective (or Laplace radius, i?off, which is 
given by lni?eff = / dqln(r). Each of P2 to P5 is con- 
sistent with having a universal non-zero asymptote, and 
moreover they are all compatible with a single common 
correction to scaling exponent 0.33, see Table I. Fig. ^ 
collects the resulting finite size scaling plots assuming 
this exponent. Together with the relative penetration 
depth results, this presents strong evidence for universal 
asymptotic geometry for DLA clusters, and a universal 
leading correction to scaling exponent v = 0.33. 

TABLE I. Best fit scaling exponents for P2 to P5. We 
have also measured Ps to Pio. These yield somewhat larger 
apparent exponents with large statistical errors. 
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FIG. 5. Finite size scaling plots for P2 to P5 for A = 1 
and A = 0.01. AH the multipole powers exhibit the same 
correction to scaling exponent. 



In all of the measurements discussed above, the cluster 
center used was the "center of charge" , natural to a snap- 
shot of the growth. In the following Section, however, 
we will require to compare data at different cluster sizes 
where it becomes natural to use a fixed center, namely the 
cluster "seed". Accordingly we have also measured the 
finite size scaling of various lengths with the seed as fixed 
origin, and in all cases using direct ensemble averages and 
for clusters with no noise reduction {A = 1). Using the 
seed as center also naturally leads to the measurement of 
penetration depth as the rms spread of deposition radius 
about its ensemble average: 



Co = Vir') - {r? 



(4) 



rather than computing the variance cluster-by-cluster be- 
fore averaging, i.e. 
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FIG. 6. Correction to scaling plots of various quantities 
of dimension length (without noise reduction, A — 1). The 
largest correction is obtained by the penetration depth ^. The 
inset magnifies the y axis around 1. 

TABLE II. Coefficients of correction to scaling fits of form 
Eq. (|), with D = 1.711 and u = 0.33. The various lengths 
are radius of deposition Rdcp = (r), seed to center of charge 



distance Rc = y {\ J dqr\^), effective radius -Rcff = exp(lnr), 
gyration radius i?gyr = \/ jf'^^i-ii^^) N' , and ensemble pen- 



etration depth ^0 = \/ (^^) — (f)'^, where the averages are over 
the ensemble of clusters at fixed A'^. 





Rdcp 


Rc 


RcS 


^gyr 


Co 


R 


0.733 


0.027 


0.726 


0.501 


0.091 


R 


-0.04 


15 


-0.14 


0.12 


6.9 



\ 



dq r"^ — [ / dqr 



(5) 



Fig. 1^ shows the fits of the following form: 

R{N) ^ RN^/°{1 + RN-") , (6) 

and the coefficients are collected in Table ||. 

It is worth noting that the effect of changing center 
is negligible except perhaps for the penetration depth, 
where for ^ using the center of of charge (as for Figs. ^ 
and H), the coefficients are | — 0.089 and ^ — 5.8 as 
opposed to those shown for in Table |l|. 

IV. MULTISCALING 

Now we consider the anomalous scaling claim of multi- 
scaling, when the aggregate has a fractal dimension which 
depends upon distance from the seed as a fraction of clus- 
ter radius. It was proposed in Ref. jsj that the ensemble 
average of the density of particles gN{r) of an A^-particle 
cluster at distance r away from the seed obeys 



gN^xRg 



A{x)R-J+^(-) 



(7) 



are 



where the dimension D{x) is function of a; = r/i?, 
and the size N and (average) radius of gyration -Rgyr 
of course mutually dependent. Using the above formula 
at fixed X, one can extract the dimension D{x) by the 
scaling with i?gyr : 



d + D{x) 



g 



(91ni?, 



yr 



(8) 



R, 



gyr 



d\ngN{xRgy,{N)) 



dRgy,/dN 
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Simple fractal scaling would require D{x) — D inde- 
pendent of X, but the dimension measured this way in 
Ref. Ql using medium size clusters (N = 10^ . . . 10^) was 
observed to be a non-trivial function (see Fig. ||). Oth- 
ers partly confirmed that findings, although with mixed 
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results 

Now we will repeat the same measurement procedure 
but instead of direct simulation we use the correction to 
scaling results of the previous Section, within a scaling 
function assumption (see below). This turns out to agree 
quantitatively with the earlier published D{x) data, but 
implies that the ultimate behaviour is simple fractal scal- 
ing with D{x) D for all x. 

Consider the distribution of r, the distance of attaching 
particles from the seed: as we have seen, this has mean 
i?dop and variance ^o- Now we assume that the shape of 
the probability density function is independent of N: 



1 



?o(iV) 



r - Rd,p{N) 
UN) 



(9) 



where ft, is a normalized probability density with zero 
mean and unit variance. After replacing the sum over 
particles with an integral, for the particle density we get 



9Nir) 



N 



dN' 
UN') 



h 



r - Ri,p{N') 
UN') 



(10) 



A similar formula has been suggested in (l3| . Given that 
we have already studied Rdep{N), Rgyr{N) and S,o{N), 
the only outstanding quantity to be found is the scaling 
function h, which we find to be very close to the standard 
normal distribution, see Fig. |^. 

Fig. ^ shows how well D{x) derived from our finite size 
scaling results plus a normal distribution for h agrees 
with the raw data of Ref. 1^. Also shown is what 
our results imply for the behaviour at larger N, and as 
N ^ CO with i?dep, ^gyr and approaching pure scal- 
ing, D{x) D. Thus we conclude that all the apparent 
reported x-dependence oi D{x) arises from corrections to 
scaling, and indeed almost all the effect cames from the 
relatively large corrections to scaling in ^o- Our new in- 
terpertation of this data also resolves a previously noted 
paradox namely that D(x) increasing with x cannot 
be asymptotic scaling as it would imply some decrease of 
(7Ar(r) with increasing N at fixed r. 
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V. SIZE FLUCTUATIONS AND FIXED POINT 



measured 
std normal 




FIG. 7. The scaling function h. The measured data at 
A'' = 10* (continuous hne) is compared to standard normal 
distribution (smooth dotted line). For the measurement, the 
histogram bin width was Au — 0.01. 



We now present an interpretation of the leading correc- 
tion to scaling, based on new observations from our data 
and building on earlier work The amplitude of the 
leading correction to scaling crosses zero at a common 
value of noise reduction Af « 0.01, for all of the plots 
in Figs. II and |[ This suggests that the noise reduction 
and the correction to scaling are fundamentally related, 
which can be understood by using the renormalization 
theory of noise reduction of Barker and Ball . In this 
view, the cluster is approximated as being at its large N 
fixed point if one unit of growth acts as a coarse graining 
of DLA on finer length scales. This seems to occur if we 
grow with input noise near Af. This is equivalent to have 
6N/N — ^/Af for relative fluctuation in the number of 
particles to advance the growth locally by one particle 
diameter. 

We can also view this in terms of a fixed point for the 
noise output of the growth, VAmt = SN/N, in terms of 
the relative fluctuation in the number of particles to span 
a fixed radius. Fig. |9| shows our data for the ensemble 
spread of extremal cluster radius. Since this spread is 
small, we can infer: 
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FIG. 8. Comparison of "multiscaling dimensions" from 
Ref. 1^ and the finite size scaling prediction discussed in the 
text. The finite size scaling prediction implies that as A'^ ^ oo, 
D{x) D for all x, and the predicted approach at A'^ = 10^ 
is shown. The only inputs to the finite size scaling curves 
are Rdm, Rgyr and using Eq. (^) with parameters from 
Table [y, plus a Gaussian model for the scaling function (see 
Fig. |). 



SN 



D 



0.060 ±0.005 



(11) 



N 



from our extrapolated value. Thus we find an asymptotic 
renormalised noise A* = 0.0036±0.0006. This is in quali- 
tative agreement with our observed value of Af. Further- 
more, Fig. |l^ shows how well this vindicates Barker and 
Ball's earlier estimates of the fixed point, using our value 
of v to extrapolate from their finite size calculations. By 
contrast, the more recent work of Cafiero et al |14| using 
a very small scale renormalisation scheme disagrees by 
two orders of magnitude. 

Our interpretation is thus that the renormalized noise 
is the slow variable which dominates convergence of other 




FIG. 9. Ensemble spread of extremal cluster radius, which 
tends towards the value 0.035 ± 0.003. 
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FIG. 10. Estimate of the fixed point value of A {A*), com- 
bined with previous estimates from Barker and Ball (the mid- 
dle three points) and Cafiero (the rightmost point). The 
Barker and Ball data is in good agreement with our results, 
but the Cafiero data disagrees. 

quantities to scaling. Our observed input noise value 
of v4/ ~ 10~^ (for the leading correction to scaling to 
vanish) and the extrapolated fixed point output noise 
A* are equal within a factor of order unity, showing the 
consistency of the picture. 

We can take this interpretation a step further to infer 
that the dominant fluctuations of i?cxt determining the 
noise reduction are fluctuations in cluster shape rather 
than overall cluster radius. The basis for this is that 
the logarithmic average radius, RcS, has much smaller 
spread, asymptotically SRcs/ReS = 0.012 ± 0.001 com- 
pared to SRcxt/Rcxt = 0.035 ± 0.003. Since i?cff is an 
average which emphasizes typical size, the larger fluc- 
tuations in i?oxt which gave us A* must be attributed 
to shape. (However, we showed in that R^s has the 
same crossover exponent, as the other quantities dis- 
cussed here.) In this sense DLA clusters are fundamen- 
tally stochastic objects with a distribution of shape. 

VI. SUMMARY 

We believe our work opens the way to a definitive view 
of DLA in two dimensions, and the extension of this work 
to three dimensions is in hand. The identification of 
'DLA fixed point behaviour' is now reasonable, as we 
have shown the sort of universal limiting amplitudes and 
correction to scaling exponents associated with such ter- 
minology. 

Some main areas are outstanding. First, the renor- 
malized noise, SN/N = ^/ A* , is not of order unity, as we 
might expect a priori, and as has been suggested @. 
We do not understand the origin of this small number, 
and tracing its origin is a central remaining challenge in 
understanding DLA. Another puzzle which we hope to 
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FIG. 11. Cluster a) grown with noise reduction A — 0.03 is 
compared with clusters grown without noise reduction, A — 1: 
b) has same gyration radius, c) has same number of particles 
and d) has same correction to scaling properties (e.g. rela- 
tive penetration depth). The numbers shown are the particle 
number and gyration radius. 

address in a later paper is why the fractal dimension is 
comparatively insensitive to the convergence of the renor- 
malized noise. 

Also, we need to understand the full scaling of the 
probability distribution for growth in DLA, correspond- 
ing to the harmonic measure of the perimeter. To 
this end the more expensive cluster growth methods of 
Hastings and Levitov |l0| (HL) are likely to come into 
their own as they yield the harmonic measure directly. 
Stepanov and Levitov ^ have already shown some re- 
sults for HL clusters grown with shallow bumps, corre- 
sponding rather closely to our noise reduction technique. 

The richer, simpler area to explore is the response to 
anisotropy and its sensitivity to noise. Small DLA clus- 
ters appear robust to the intrinsic bias of growing on a 
square lattice, whereas large clusters (and equivalently 
noise reduced ones) are driven to grow a four fingered 
dendrite. The first requirement is a systematic analysis of 
how this is a relevant perturbation of the isotropic DLA 
fixed point. Secondly, we might ask whether the anoma- 
lous response for small simple DLA clusters is domi- 
nated by some other hitherto unsuspected fixed point 
with much larger noise level. There is another rather 
neglected nearby fixed point, that of spherical growth, 
which becomes more pertinent at high noise reduction - 
where it takes longer to exhibit its instability. We suggest 
the influence of this flxed point may be responsible for 
shifting the observed Af somewhat above A*, and this 
should be relatively amenable to analytic theory. 



6 



ACKNOWLEDGMENTS 



NEB would like to thank BP Amoco and EPSRC for 

the support of a CASE award during this research. ES is 
supported by the Dutch FOM Foundation and the EU 
Marie Curie Fellowship. We thank Paul Meakin and 
Thomas Rage for sending us a computer code which we 
used for some of the results presented here. 



[1] T. A. Witten and L. M. Sander, Phys. Rev. Lett. 47, 

1400 (1981). 

[2] E. Somfai, L. M. Sander, and R. C. Ball, Phys. Rev. Lett. 

83, 5523 (1999). 
[3] M. Plischke and Z. Racz, Phys. Rev. Lett. 53, 415 (1984). 

[4] B. Davidovitch, H. G. E. Hcntschcl, Z. Olami, L Procac- 
cia, L. M. Sander, and E. Somfai, Phys. Rev. E 59, 1368 
(1999). 

[5] A. Coniglio and M. Zannetti, Physica A 163, 325 (1990). 
[6] C. Amitrano, A. Coniglio, P. Meakin, and M. Zannetti, 

Phys. Rev. B 44, 4974 (1991). 
[7] C. Tang, Phys. Rev. A 31, 1977 (1985). 
[8] P. W. Barker and R. C. Ball, Phys. Rev. A 42, 6289 

(1990). 

[9] M. G. Stepanov and L. S. Levitov, Phys. Rev. E 63, 
061102 (2001). 

[10] M. B. Hastings and L. S. Levitov, Physica D 116, 244 
(1998). 

[11] P. Ossadnik, Physica A 176, 454 (1991). 
[12] P. Ossadnik, Physica A 195, 319 (1993). 
[13] J. Lee, S. Schwarzer, A. Coniglio, H. E. Stanley, 

Phys. Rev. E 48, 1305 (1993). 
[14] R. Cafiero, L. Pietronero, and A. Vespignani, Phys. Rev. 

Lett. 70, 3939 (1993). 



7 



